MRI APPARATUS 



BACKGROUND OF THE INVENTION 

The present invention relates to an MR (magnetic resonance) image 
producing method and an MRI (magnetic resonance imaging) apparatus, and 
more particularly to an MR image producing method and MRI apparatus that 
improve on rendering capability for a blood vessel. 

A conventional MRI apparatus comprises MR signal acquiring means for 
acquiring MR signals, window-processing means for window-processing the MR 
signals using a window function that has a value of "one" from the center of a 
k-space to a position near the periphery of the k-space and has a decreasing value 
as it goes closer to the periphery, and Fourier-transformation processing means 
for conducting Fourier-transformation processing on the window-processed MR 
signals to produce an MR image. 

The window processing is conducted to concentrically suppress a high 
frequency portion of MR signals to thereby prevent truncation artifacts or 
anisotropic noise texture due to signal acquisition confined to a limited 
rectangular region in a k-space by the MRI apparatus. 

Related conventional techniques are disclosed in Japanese Patent 
Application Laid Open Nos. H4-53531 and H6-121781. 

In the conventional MRI apparatus, the same window processing is 
conducted whether the image to be produced is a blood flow image or not. 

In other words, the conventional window processing is not optimal when 
the image to be produced is a blood flow image, and it cannot improve rendering 
capability for a blood vessel. 

SUMMARY OF THE INVENTION 

It is therefore an object of the present invention to provide an MR image 
producing method and MRI apparatus that improve on rendering capability for a 
blood vessel by optimizing window processing when an image to be produced is 
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a blood flow image. 

In a first aspect, the present invention provides an MR image producing 
method characterized in comprising: window-processing MR signals using a 
window function that has a value less than one at a center and in its proximate 
5 region in a k-space and on a periphery and in its proximate region in the k-space, 
and, between the regions in which the window function has a value less than one, 
has a value larger than that in the regions in which the window function has a 
value less than one; and applying Fourier-transformation processing to the 
window-processed MR signals to obtain an MR image. 

10 In this configuration, the proximate region of the center of the k-space is a 

range of about five — twenty data points from the center of the k-space. The 
proximate region of the periphery of the k-space is a range of about five — 
twenty data points from the periphery of the k-space. 

According to the MR image producing method of the first aspect, since a 

15 window function that has a value less than one at a center and in its proximate 
region in the k-space is employed, MR signals near the center of the k-space are 
suppressed. MR signals of a tissue portion are narrowly distributed near the 
center of the k-space, while MR signals of a blood flow portion are broadly 
distributed in a high frequency region, as well as near the center. Thus, the MR 

20 signals of the tissue portion are strongly suppressed, while the MR signals of the 
blood flow portion are relatively weakly suppressed. Therefore, rendering 
capability for a blood vessel is relatively improved. 

Moreover, since the window function has a value less than one on a 
periphery and in its proximate region in the k-space, a high frequency portion of 

25 MR signals can be concentrically suppressed as in the conventional technique. 

In a second aspect, the present invention provides an MR image producing 
method characterized in comprising: window-processing MR signals using a 
window function that has a value less than one at a center of a k-space, first 
increases to a value C equal to or more than one as it goes farther from the center, 

30 remains at C for a certain duration, then passes to one, and decreases to a value 
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less than one as it goes from near a periphery to the periphery of the k-space; and 
applying Fourier-transformation processing to the window-processed MR signals 
to obtain an MR image. 

In this configuration, the region in which the window function has a value 
5 increasing to a value C from the center of the k-space is a range of about three — 
fifteen data points from the center of the k-space. The region in which the 
window function remains at C for a certain duration is a range of about twenty 
— fifty data points. The region in which the window function passes from C to 
one is a range of about three — ten data points. The region in which the 

10 window function has a value less than one is a range of about five — twenty data 
points from the periphery of the k-space. 

According to the MR image producing method of the second aspect, since 
the window function has a value less than one at a center and in its proximate 
region in the k-space, MR signals near the center of the k-space are suppressed. 

15 MR signals of a tissue portion are narrowly distributed near the center of the 
k-space, while MR signals of a blood flow portion are broadly distributed in a 
high frequency region, as well as near the center. Thus, the MR signals of the 
tissue portion are strongly suppressed, while the MR signals of the blood flow 
portion are relatively weakly suppressed. Next, in the region in which "the 

20 window function remains at C for a certain duration", a zero-th order peak 
portion (a crest having a maximum value at the center) in the MR signals of the 
blood flow portion is preserved or amplified. Next, in the region in which "the 
window function passes to one", a first- or higher order peak portion (a crest 
having a maximum value at a position except the center) in the MR signals of the 

25 blood flow portion is preserved. Thus, rendering capability for a blood vessel is 
relatively improved. 

Moreover, since the window function has a value less than one on a 
periphery and in its proximate region in the k-space, a high frequency portion of 
MR signals can be concentrically suppressed as in the conventional technique. 

30 In a third aspect, the present invention provides the MR image producing 
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method having the aforementioned configuration, characterized in that: the 
window function is a function using a Gaussian function in the region in which 
the window function increases to C. 

According to the MR image producing method of the third aspect, a 
5 Gaussian function exp{- | k | 2 /a 2 } can be used to smoothly increase the value 
from a value less than one to a value C. 

In a fourth aspect, the present invention provides the MR image producing 
method having the aforementioned configuration, characterized in that: the 
window function is a function using a Fermi-Dirac function in the region in 

10 which the window function decreases to a value less than one. 

According to the MR image producing method of the fourth aspect, a 
Fermi-Dirac function 1 / (1 + exp{( | k | - R) / b}) can be used to smoothly reduce 
the value from one to a value less than one. 

In a fifth aspect, the present invention provides an MR image producing 

15 method characterized in comprising: window-processing MR signals using a 
window function that has a value less than one at a center of a k-space, first 
increases to one as it goes farther from the center, remains at one for a certain 
duration, and decreases to a value less than one as it goes from near a periphery 
to the periphery of the k-space; and applying Fourier-transformation processing 

20 to the window-processed MR signals to obtain an MR image. 

According to the MR image producing method of the fifth aspect, since the 
window function has a value less than one at a center and in its proximate region 
in the k-space, MR signals near the center of the k-space are suppressed. MR 
signals of a tissue portion are narrowly distributed near the center of the k-space, 

25 while MR signals of a blood flow portion are broadly distributed in a high 
frequency region, as well as near the center. Thus, the MR signals of the tissue 
portion are strongly suppressed, while the MR signals of the blood flow portion 
are relatively weakly suppressed. Next, in the region in which "the window 
function remains at one for a certain duration/' MR signals of the blood flow 

30 portion are preserved. Thus, rendering capability for a blood vessel is relatively 
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improved. 

Moreover, since the window function has a value less than one on a 
periphery and in its proximate region in the k-space, a high frequency portion of 
MR signals can be concentrically suppressed as in the conventional technique. 
5 In a sixth aspect, the present invention provides the MR image producing 

method having the aforementioned configuration, characterized in that: the 
window function is a function using a Gaussian function in the region in which 
the window function increases to one. 

According to the MR image producing method of the sixth aspect, a 
10 Gaussian function exp{- | k 1 2 / a2 } can be used to smoothly increase the value 
from a value less than one to one. 

In a seventh aspect, the present invention provides the MR image 
producing method having the aforementioned configuration, characterized in 
that: the window function is a function using a Fermi-Dirac function in the region 
15 in which the window function decreases to a value less than one. 

According to the MR image producing method of the seventh aspect, a 
Fermi-Dirac function 1 / (1 + exp{( | k | - R) / b}) can be used to smoothly reduce 
the value from one to a value less than one. 

In an eighth aspect, the present invention provides an MR image producing 
20 method characterized in comprising: producing MR images by the MR image 
producing method having the aforementioned configuration for a plurality of 
sequential slices; generating three-dimensional data from the MR images; and 
conducting MIP processing on the three-dimensional data to produce a 
projection image. 

25 In the MR image producing method of the eighth aspect, rendering 

capability for a blood vessel can be improved for an angiographic image. 

In a ninth aspect, the present invention provides an MRI apparatus 
characterized in comprising: MR signal acquiring means for acquiring MR 
signals; window-processing means for window-processing the MR signals using 

30 a window function that has a value less than one at a center and in its proximate 
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region in a k-space and on a periphery and in its proximate region in the k-space, 
and, between the regions in which the window function has a value less than one, 
has a value larger than that in the regions in which the window function has a 
value less than one; and Fourier-transformation processing means for applying 
5 Fourier-transformation processing to the window-processed MR signals to obtain 
an MR image. 

According to the MRI apparatus of the ninth aspect, the MR image 
producing method of the first aspect can be suitably implemented. 

In a tenth aspect, the present invention provides an MRI apparatus 
10 characterized in comprising: MR signal acquiring means for acquiring MR 
signals; window-processing means for window-processing the MR signals using 
a window function that has a value less than one at a center of a k-space, first 
increases to a value C equal to or more than one as it goes farther from the center, 
remains at C for a certain duration, then passes to one, and decreases to a value 
1 5 less than one as it goes from near a periphery to the periphery of the k-space; and 
Fourier-transformation processing means for applying Fourier-transformation 
processing to the window-processed MR signals to obtain an MR image. 

According to the MRI apparatus of the tenth aspect, the MR image 
producing method of the second aspect can be suitably implemented. 
20 In an eleventh aspect, the present invention provides the MRI apparatus 

having the aforementioned configuration, characterized in that: the window 
function is a function using a Gaussian function in the region in which the 
window function increases to C. 

According to the MRI apparatus of the eleventh aspect, the MR image 
25 producing method of the third aspect can be suitably implemented. 

In a twelfth aspect, the present invention provides the MRI apparatus 
having the aforementioned configuration, characterized in that: the window 
function is a function using a Fermi-Dirac function in the region in which the 
window function decreases to a value less than one. 
30 According to the MRI apparatus of the twelfth aspect, the MR image 
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producing method of the fourth aspect can be suitably implemented. 

In a thirteenth aspect, the present invention provides an MRI apparatus 
characterized in comprising: MR signal acquiring means for acquiring MR 
signals; window-processing means for window-processing the MR signals using 
5 a window function that has a value less than one at a center of a k-space, first 
increases to one as it goes farther from the center, remains at one for a certain 
duration, and decreases to a value less than one as it goes from near a periphery 
to the periphery of the k-space; and Fourier-transformation processing means for 
applying Fourier-transformation processing to the window-processed MR signals 
10 to obtain an MR image. 

According to the MRI apparatus of the thirteenth aspect, the MR image 
producing method of the fifth aspect can be suitably implemented. 

In a fourteenth aspect, the present invention provides the MRI apparatus 
having the aforementioned configuration, characterized in that: the window 
15 function is a function using a Gaussian function in the region in which the 
window function increases to one. 

According to the MRI apparatus of the fourteenth aspect, the MR image 
producing method of the sixth aspect can be suitably implemented. 

In a fifteenth aspect, the present invention provides the MRI apparatus 
20 having the aforementioned configuration, characterized in that: the window 
function is a function using a Fermi-Dirac function in the region in which the 
window function decreases to a value less than one. 

According to the MRI apparatus of the fifteenth aspect, the MR image 
producing method of the seventh aspect can be suitably implemented. 
25 In a sixteenth aspect, the present invention provides the MRI apparatus 

characterized in comprising: three-dimensional data generating means for 
generating three-dimensional data from MR images produced for a plurality of 
sequential slices; and MlP-processing means for conducting MIP processing on 
the three-dimensional data to produce a projection image. 
30 According to the MRI apparatus of the sixteenth aspect, the MR image 
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producing method of the eighth aspect can be suitably implemented. 

According to the MR image producing method and MRI apparatus of the 
present invention, rendering capability for a blood vessel is improved. 
Moreover, noise in a high frequency region is suppressed to improve CNR 
5 (carrier to noise ratio). 

Further objects and advantages of the present invention will be apparent 
from the following description of the preferred embodiments of the invention as 
illustrated in the accompanying drawings. 

1 0 BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 is a block diagram showing an MRI apparatus in accordance with 
a first embodiment. 

Figure 2 is a flow chart showing MR image producing processing in 
accordance with the first embodiment. 
15 Figure 3 is a graph showing the MR signal intensity on a center axis of a 

slice in an assumed model. 

Figure 4 is a graph showing the MR signal intensity on a center axis of a 
k-space in the assumed model. 

Figure 5 is a graph showing a value of a window function on a center axis 
20 of the k-space in accordance with the first embodiment. 

Figure 6 is a graph showing the MR signal intensity on a center axis of the 
k-space after window processing in accordance with the first embodiment. 

Figure 7 is a graph showing the MR signal intensity on a center axis of a 
slice after FFT processing in accordance with the first embodiment. 
25 Figure 8 is a graph showing a value of a window function on a center axis 

of a k-space in accordance with a second embodiment. 

Figure 9 is a flow chart showing three-dimensional MR image processing in 
accordance with a third embodiment. 



30 DETAILED DESCRIPTION OF THE INVENTION 
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The present invention will now be described in more detail with reference 
to embodiments shown in the accompanying drawings. 
- First Embodiment - 

Figure 1 is a block diagram showing an MRI apparatus in accordance with 
5 a first embodiment. 

In the MRI apparatus 100, a magnet assembly 1 has a bore (cavity portion) 
for inserting therein a subject, and is provided, surrounding the bore, with a 
gradient coil (which comprises X-axis, Y-axis and Z-axis coils, and the 
combination thereof determines slice, warp and read axes) 1G for generating 
10 gradient magnetic fields, a transmit coil IT for applying RF pulses for exciting 
spins of atomic nuclei within the subject, a receive coil 1R for detecting NMR 
signals from the subject, and a static magnetic field power supply 2 and static 
magnetic field coil 1C for generating a static magnetic field. 

It should be noted that permanent magnets may be employed in place of 
15 the static magnetic field power supply 2 and static magnetic field coil 1C 
(superconductive coil). 

The gradient coil 1G is connected to a gradient coil driving circuit 3. The 
transmit coil IT is connected to an RF power amplifier 4. The receive coil 1R is 
connected to a preamplifier 5. 
20 A sequence memory circuit 8 operates the gradient coil driving circuit 3 

based on a stored pulse sequence in response to instructions from a computer 7 
to thereby generate gradient magnetic fields from the gradient coil 1G. The 
sequence memory circuit 8 also operates a gate modulation circuit 9 to modulate 
high frequency output signals from an RF oscillation circuit 10 into pulsed 
25 signals of predefined timing and envelope. The pulsed signals are applied to 
the RF power amplifier 4 as excitation pulses, power-amplified in the RF power 
amplifier 4, and then applied to the transmit coil IT in the magnet assembly 1 to 
transmit RF pulses. 

The preamplifier 5 amplifies NMR signals from the subject detected at the 
30 receive coil 1R in the magnet assembly 1, and inputs the signals to a phase 
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detector 12. The phase detector 12 phase-detects the NMR signals from the 
preamplifier 5 employing the output from the RF oscillation circuit 10 as a 
reference signal, and supplies the phase-detected signals to an A/D converter 11. 
The A/D converter 11 converts the phase-detected analog signals into MR data in 
the form of digital signals, and inputs them to the computer 7. 

The computer 7 reads the MR data from the A/D converter 11, and 
performs image reconstruction processing to produce an MR image. The 
computer 7 is also responsible for overall control such as receiving information 
supplied from an operator console 13. 

A display device 6 displays the MR image. 

Figure 2 is a flow chart showing MR image producing processing by the 
MRI apparatus 100. 

At Step PI, window processing is applied to MR data collected from one 

slice. 

At that time, a model for the slice is assumed such that a circular blood 
vessel lies in the center of a 256 x 256 FOV (field of view) and brain matters 
surround the blood vessel in a donut-like shape. Moreover, for the MR signal 
intensity on a slice center axis, the model is assumed to have the maximum signal 
intensity of the brain matters G and that of the blood vessel V equal to each other, 
as shown in Figure 3. The horizontal axis in Figure 3 represents the index of a 
pixel point. 

The MR signal intensity over a whole slice plane is a curved surface formed 
by axially rotating the curve of Figure 3 around an axis of rotation at a position O 
in Figure 3. 

For this model, a distribution of frequency components of MR signals on 
the center axis of the k-space is shown in Figure 4. The horizontal axis in Figure 
4 represents the index of a data point in the k-space. From the viewpoint of 
frequency, a position O on the horizontal axis corresponds to a direct-current 
position (the center of the k-space), and the frequency becomes higher as it 
approaches the periphery. 
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As shown in Figure 4, the frequency component g of the brain matters G 
contained in the MR data has a very strong peak at the center, and is narrowly 
distributed near the center. On the other hand, the frequency component v of 
the blood vessel V (the frequency component of blood flow) contained in the MR 
5 data has a relatively weak peak at the center, and is broadly distributed in a high 
frequency region, as well as near the center. 

The distribution of frequency components of MR signals over a whole 
k-space is a curved surface formed by axially rotating the curve of Figure 4 
around an axis of rotation at the center in Figure 4. 
10 Figure 5 shows a window function f(k) employed for MR signals on the 

center axis of the k-space. 

The window function f(k) has a value 0.5 at the center of the k-space, first 
increases to a value 1.2 as it goes farther from the center, remains at the value 1.2 
fox a certain duration, then passes to one, and decreases to a value 0.5 as it goes 
15 from near a periphery to the periphery of the k-space. 

In the region in which the value increases from 0.5 to 1.2, the window 
function f (k) is represented by: 

f(k) = C(l-A.exp{-|k|Va2}), 
where a, A and C are constants, and k is the distance from the center. In Figure 
20 5,C = 1.2, A = 0.6, and a = 5. 

In the region in which the value decreases from one to 0.5, the window 
function f (k) is represented by: 

f(k) = l/(l + exp{(|k| -R)/b}), 
where b and R are constants, and k is the distance from the center. In Figure 5, R 
25 = 128, and b = 3. 

The window function f (k) over a whole k-space is a curved surface formed 
by axially rotating the curve of Figure 5 around an axis of rotation at the center in 
Figure 5. 

Figure 6 shows a distribution of frequency components in the 
30 window-processed MR signals on a center axis of the k-space. 
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The frequency component g of the brain matters G is strongly suppressed 
to about a half of the original intensity. 

On the other hand, a zero-th order peak portion of the frequency 
component v of the blood vessel is suppressed at the center, while it is amplified 
near the center. A first- and higher order peak portions are preserved at the 
original intensity. 

The frequency components of MR signals over a whole k-space are a 
curved surface formed by axially rotating the curve of Figure 6 around an axis of 
rotation at the center in Figure 6. 

Returning to Figure 2, at Step P2, the window-processed MR data is 
processed by two-dimensional FFT (fast Fourier transformation). 

Figure 7 shows the MR signal intensity on the center axis of a slice after the 
FFT processing. 

The signal of the brain matters G is strongly suppressed, while the signal of 
the blood vessel V remains near the original intensity. 

The MR signal intensity over a whole slice is a curved surface formed by 
axially rotating the curve of Figure 7 around an axis of rotation at a position O in 
Figure 7. 

As a result, the rendering capability for the blood vessel V can be relatively 
improved. 

Moreover, the window function f(k) has a value less than one on the 
periphery and in its proximate region in the k-space. Thus, a high frequency 
portion of MR signals can be concentrically suppressed as in the conventional 
technique. 

- Second Embodiment - 

In a second embodiment, a window function as shown in Figure 8 is used 
to achieve window processing. 

The window function f (k) in Figure 8 has a value 0.5 at the center of the 
k-space, increases to a value one as it goes farther from the center, remains at the 
value one for a certain duration, and decreases to a value 0.5 as it goes from near 
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the periphery of the k-space to the periphery. 

In the region in which the value increases from 0.5 to one, the window 
function f (k) is represented by: 

f(k) = l-A.exp{-|k|2/ a 2}, 
where a and A are constants, and k is the distance from the center. 

In the region in which the value decreases from one to 0.5, the window 
function f (k) is represented by: 

f(k)-l/(l+exp{(|k| -R)/b}), 
where b and R are constants, and Zeis a distance from the center. 

The window function f (k) over a whole k-space is a curved surface formed 
by axially rotating the curve of Figure 8 around an axis of rotation at the center in 
Figure 8. 

It should be noted that the window function would be identical to the 
conventional one when it has a value of one at the center and in its proximate 
region, as indicated by broken line in Figure 8. 
- Third Embodiment - 

In a third embodiment, a three-dimensional MR image is produced with 
improved rendering capability for the blood vessel V. 

Figure 9 is a flow chart showing three-dimensional MR image producing 
processing in accordance with the third embodiment. 

At Step Ql, an MR image is produced according to one of the 
aforementioned embodiments, and the production is repeated for a plurality of 
slices sequentially lining up in a thickness direction. 

At Step Q2, three-dimensional data is generated from the MR images of the 
plurality of slices sequentially lining up in the thickness direction. 

At Step Q3, MIP (maximum intensity projection) processing is conducted 
on the three-dimensional data to produce a three-dimensional MR image. 

Since the MR images of the slices are those with improved rendering 
capability for the blood vessel V, the three-dimensional MR image is also an MR 
image with improved rendering capability for the blood vessel V, and the present 
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embodiment is suitable for producing an angiographic image. 
- Other Embodiments - 

The window function may be generated from a combination of functions 
other than the aforementioned functions for use in the window processing. 

Many widely different embodiments of the invention may be configured 
without departing from the spirit and the scope of the present invention. It 
should be understood that the present invention is not limited to the specific 
embodiments described in the specification, except as defined in the appended 
claims. 



